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Many recent works on large A'^ holographic QCD in the planar limit have not considered UV 
completions, restricting exclusively towards analyzing the IR physics. Due to this, the UV problems 
like Landau poles and divergences of Wilson loops including instabilities at high temperatures have 
not been addressed. In some of our recent papers, we have discussed a possible UV completion, 
which is conformal in the UV and confining in the far IR, that avoids the Landau poles and the 
Wilson loop divergences. In this paper we give a field theory realization of this including the 
complete RG flow. We extend our UV complete model to study scenarios both above and below 
the deconfinement temperature and argue how phase transition in our model should be understood. 
Interestingly, because of the UV completion, subtle issues like instability due to negative specific 
heat do not appear. We also briefly elucidate the advantages that our model may have over other 
models studying large A'^ thermal QCD. 



1. INTRODUCTION 

The gauge/gravity duality has so far proved to be a 
powerful technique to solve many strong coupling prob- 
lems of large N gauge theories, and especially large N 
QCD, in the planar limit. The application of a gravity 
dual to understand strongly coupled gauge theory was, 
in retrospect, the next best thing to do. A simple way to 
see this would be to consider a particular gauge-theory 
defined on a 3 -I- 1 dimensional slice at a certain energy 
scale A. Now imagine we stack up all the slices together, 
described at different energy scales, along an orthogonal 
direction (call it the "radial" direction r). This way we 
will get a five dimensional space that captures the full 
dynamics of a given gauge theory from the Ultra- Violet 
(UV), i.e large r, to the Infra- Red (IR), i.e small r. The 
"radial" direction would then obviously be the direction 
along which the energy would change, i.e the direction of 
the Renormalisation Group (RG) flow. For a Conformal 
Field Theory (CFT), the theory does not change along 
the radial direction^ and therefore could as well be de- 
fined at the boundary of the five-dimensional space. The 
scale invariance of the underlying gauge theory will re- 
strict the geometry of the five-dimensional space to the 
Anti-deSitter (AdS) space 1], although it would be in- 
teresting to argue that this is the unique choice^. 

However, for gauge theories with inherent RG flows. 
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^ Assuming the usual behavior of the irrelevant operators. 

^ Furthermore, a Feynman diagram for any interaction between 
point-like particles, when stacked up as above, would look like 
an interaction between extended objects, i.e strings! This is ba- 
sically the essence of using string (or gravity) duals to study 
gauge theories. It will be informative to make this more precise. 



the situation will be different and it would be instructive 
to study the theories at various r (although we could also 
restrict ourselves to the boundary again). The example 
that we are interested in is large N QCD, which we expect 
to be asymptotically conformal^ in the UV and confining 
in the far IR. Specific geometries that do the jobs for 
both zero and non-zero temperatures were presented in 
Q although the details of the gauge theories were not 
presented there. In this paper we will fill up some of the 
gaps left in d, @ and argue why we believe our choice of 
the gravity dual is better suited to study large N thermal 
QCD (see also IT] for another model that studies UV 
complete large N thermal QCD from a bottom-up five- 
dimensional point of view). 



2. THE FIELD THEORY FROM THE GRAVITY 
DUAL 

The gravity dual of a large N thermal QCD above the 
deconfinement temperature, described using only a fla- 
vored Klebanov-Strassler geometry with a black-hole 
has few ultra-violet (UV) problems. For example, there 
are Landau poles coming from the flavor branes, and 
the Wilson loops are generically UV divergent Q. All 
these issues could be resolved if we properly augment the 
Klebanov-Strassler geometry, which we will henceforth 
call as the Ouyang-Klebanov-Strassler black-hole (OKS- 
BH) [1, @ geometry, with a suitable asymptotically 
Anti-de Sitter (AdS) space. As discussed in [6] , this aug- 
mentation can only be performed in the presence of an 



^ It is interesting that we demand conformal behavior in the UV 
and not asymptotic freedom. This is because the 'tHooft cou- 
pling A = Qyj^jN approaches a constant in the limit Qym ~^ ^ 
and N oo. This way, the theory is actually asymptotically 
free in terms of Qym conformal in terms of A. Furthermore, 
we will demand A to be very large throughout the whole RG flow 
so that the gravity dual can be restricted to its classical limit. 
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interpolating space and certain number ol anti five-brane 
sources. 

The interpolating region, which we called region 2 in 
0, can be interpreted alternatively as the deformation 
of the neighboring geometry once we attach an AdS cap 
to the OKS-BH geometry. The OKS-BH geometry is in 
the range < r < r^in (which we will call as region 
1) and the AdS cap is the range r > vq (which we will 
call as region 3). Here r/j is the horizon radius. The 
geometry in the range < r < is the deformation. 
Such deformations should be expected for all other UV 
caps advocated in la|. This construction was elaborated 
in some details in [6|. In this paper we will start with a 
gauge theory interpretation of background. 



confinement at zero temperature. _ 

The few calculations that we did in [Q] regarding (a) 
the flow of N and M colors, (b) the RG flows, (c) the 
decay of the three-forms and (d) the behavior of the dual 
gravity background all support the gauge theory interpre- 
tation that we gave above. What we haven't been able 
to demonstrate in 0, |^ is the precise Higgsing that takes 
us to the cascading picture. From the gravity side, it is 
clear how this could be interpreted. From the gauge the- 
ory side, we will provide a brief derivation below. But 
before we dwell on the details, let us see how the full 
Renormalisation Group (RG) flow would look like with 
the AdS cap. 




Smooth flow at the center 
of the RG fixet] point surface 



Energy scale 



(strongly conpletl gauge theory) 

(a) RG flow in the far (b) RG flow in the 
IR far UV 

FIG. 1: In the far IR, one may note that the cascading 
RG flow is never captured by the classical gravity theory. 
The classical supergravity description would capture only the 
smooth parts of the RG flow shown in (a) at the center of 
each slices. The vertical distances in (a) refer to the slices de- 
scribed using appropriate Seiberg dual descriptions. On the 
other hand, in (b) the RG flows all tend to go to zero at some 
UV scales. This is where the theory become conformal (all 
scales are chosen with a' = 1). Note that, only the strongly 
coupled parts of figure (b) are captured by the classical su- 
pergravity description. 



For the UV region r > tq we expect the dual gauge 
theory to be SU{N + M) x SU{N+M) with fundamental 
flavors coming from the seven-branes. This is because 
addition of M anti five-branes at the junction (i.e for 
r > ro) with gauge fluxes on its world- volume tells us that 
the number of three-branes degrees of freedom are N+M , 
where the M and N factors come from the presence of 
M five-branes anti-five-branes pairs and N D3-branes. 
Furthermore, the SU{N+M) x SU{N+M) gauge theory 
informs us that the gravity dual is approximately AdS, 
but has RG flows because of the fundamental flavors. In 
other words, the two couplings gi and g2 of each gauge 
group would be approximately the same and exhibit a 
walking RG flow. At the scale r — ro, we expect one of 
the gauge group to be Higgsed, so that we are left with 
SU{N + M) X SU{N). Now both gauge couplings flow 
at different rates and give rise to a cascade that is slowed 
down by the Nf flavors. In the end, at far IR, we expect 



2.1. Continuous RG Flow from UV to IR 

As mentioned above, the gravity dual should give us a 
RG flow that allows us to see the UV conformal behavior 
and the IR confining behavior succinctly. However, there 
is a subtlety as shown in fig 1. The cascading RG fiow in 
the far IR, where the theory goes from one Seiberg fixed 
point to another, is in fact not seen in the dual gravity 
side because it runs between weakly coupled theories. 
Thus, what we see from the dual gravity side is a smooth 
RG fiow^, as depicted at the center of each slice in fig 
1(a). 

The RG fiow in the intermediate region, identified as 
region 2, is more involved and will be discussed in details 
in [ll|. However this RG flow connects smoothly to the 
RG fiow in the AdS cap, called as region 3. The fiow in 
region 3 approaches conformality where both couplings 
run at an equal rate as shown in fig 1(b). 

The Beta functions are also easy to compute to first 
order in g^iV/ from the gravity dual. In Region 1 the two 
couplings at a scale A run in the following way: 
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(2.1) 



47r 



log(A^ + 9a 



where the RHS of both equations is evaluated at r = A 
in the gravity picture. The constant a appearing above is 
the bare resolution parameter that one may set to zero^. 
In this limit, the RG fiow is clearly the NSVZ RG fiow 



This also means that at any given scale A there are in prin- 
ciple an inflnite number of gauge theory descriptions available. 
Out of which, one of them might be the most useful description 
at that scale and is therefore captured by the classical super- 
gravity analysis at r = A. For example, in the far IR, out of 
the many available gauge theory descriptions, it is the conflning 
SU{M) gauge theory (which is naturally strongly coupled) that 
is captured by the classical supergravity solution at small r. As 
a consequence, we expect that the definition of the number of 
colors at any given scale would become a little ambiguous. 
This is however not so above the deconflnement temperature. 
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FIG. 2: A slightly unconventional way to represent the RG 
flow in our model. We get the complete RG flow by glu- 
ing the three regions altogether and using S-duality to trans- 
mute strong coupling into weak coupling. Starting from IR 
regime, once a particular coupling gets strong, a S-duality is 
performed to reverse the sign of the beta function associated 
with that coupling. This appears as the sharp edges in the 
figure above. From the UV region this can be seen in the 
following way: The coupling starts as a constant in Region 3, 
when it gets to the transition point ro = 200, it has a small 
plateau region continuing in region 2, then it flows down to 
Region 1. The RG flow continues after the transition point 
Tmin ~ 100, but siucc the rate of change is fast more sharp 
corners appear in region 1. These are the points connected to 
their S-dual values. Eventually this reaches the smallest en- 
ergy possible after which we expect linear confinement at low 
temperatures. As before, all scales are chosen with a' = 1. 



|12| . On the other hand, in region 2, where we still have 
two couplings, the RG flow is highly non-trivial. This can 
be derived from the gravity dual where we see that the 
three-form fluxes play an important role in the running 
of the couphngs [llj : 

^ = + B,] (2.2) 

gi 2 27r Jg2 

92^ 2 27r Jg2 

Finally in region 3, the scenario is somewhat simpler. 
The two couplings flow approximately at the same rate 
and the flow is governed by the Nf D7 and anti-D7 
pairs that we keep in region 3 to cancel the Landau 
poles. These seven-branes are responsible for restoring 
the SU(Nf) X SU{Nf) chiral symmetry above the de- 
confinement temperature (i.e when we insert a black- hole 
with a horizon radius 0, Q ) . The running of the cou- 
pling, which we call gvM, is now: 

^^=3YM 2^j^ (2.4) 

n— 1 

where Vn are all independent of A and whose precise 
form will be derived in [ll[ . 



As shown recently in @|, even if we demand a vanishing bare 
resolution parameter, it'll get a contribution from the horizon 
radius r/j, such that a ~ 0{rh). Of course, on the gauge theory 
side, the branes are still wrapped on vanishing cycle. 



2.2. Higgsing 

In region 3, we have a SU{N + M) x SU{N + M) 
gauge group which breaks down to SU{N + AI) x SU{N) 
by the Higgs mechanism as we enter Region 1. We will 
study this mechanism in two versions: supersymmetric 
and non-supersymmetric. Since the purpose is to break 
the gauge group, we will ignore any fundamental matter 
fields in the following discussion. 

Before moving ahead, let us see how we could jus- 
tify the Higgs mechanism from the gravity perspective. 
The brane construction that reproduces the gauge theory 
should be understood on a scale-by-scale basis, so that 
the full RG flow could be reproduced in the gravity dual. 
Generically, we expect D3s and M wrapped D5s on 
a vanishing two-cycle of the conifold. Allowing a small 
resolution factor to the other two-cycle, we can distribute 
the anti-D5 branes on the resolved sphere such that they 
wrap the same vanishing two-cycle but are distributed on 
the other sphere. Similarly the D7 and anti-D7 branes 
are also distributed^, over the resolved sphere via the 
Ouyang embedding [6]. This configuration is more intu- 
itive from the gravity dual side where the radial coordi- 
nate now becomes the scale of the theory. At a given 
scale we expect Me number of wrapped anti-D5 branes 
where = "^^^^ ^ and a » 1. 

Then it is easy to see that the resulting gauge group be- 
comes SU{N + Me) X SU{N + M). Clearly, in the limit 
e — >■ 0, we recover the conformal gauge group. There- 
fore, the anti-D5 branes appear to only affect one of the 
gauge groups in the product. In region 3, where r » tq, 
A/f w M , this tells us that the RG flow will be mostly 
due to the flavor seven-branes^. 



° The tachyons between D5 and anti D5-branes or between D7 
and anti D7-branes can be cancelled by switching on appropri- 
ate gauge fluxes on the set of anti branes. This phenomena is 
somewhat similar to the ones in [T3l . These gauge fluxes will 
create bound D3 and bound D5-branes respectively on the two 
set of brane anti-brane systems. If oriented properly, the sys- 
tem would then be almost BPS in the zero-temperature case 
when the distance between the two set of branes is large (the 
multipole forces are heavily suppressed). To stabilize this com- 
pletely, one may switch on three-form //ns i ^^RR fluxes on the 
internal space (the axio-dilaton are already switched on). These 
i?-fluxes would not only change the moding of the strings be- 
tween the branes but also stabilize the position of the branes, 
by generating perturbative and non-perturbative superpotential 
and giving masses to the scalar fields on the branes, along the 
lines of Il4l Il5l, [2^ . Alternatively for short distances, one may 
dissolve the anti-D5 branes in the D7 anti-D7 system in the way 
discussed in 0], and then stabilize the seven-brane positions. In 
either case the physics would be the same. 

The gauge group that actually appears in the far IR is SU (N + 
M) X SU{N) X U{1)^'' , where the (/(l)'s are from the massless 
sector of the anti-D5 branes. However at low energies, once we 
integrate out the Higgs masses (i.e the strings between the D5 
and the anti-D5 branes), these (7(l)'s would be decoupled. Fur- 
thermore at strong coupling, where we expect the dual gravity 
description to hold, these (/(l)'s will never appear. This in turn 
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The above construction then instructs us that the Hig- 
gsing process generating the cascade should simply be en- 
gineered by making some anti-D5 brane DOFs heavy, i.e 
by moving the anti-D5 branes away from the N D3 and 
the M wrapped D5 branes on the resolved sphere as we 
discussed above. In a supersymmetric theory, where the 
UV completion is done by a A/" = 2 theory, this process 
would mean moving the anti-D5 brane DOFs along the 
Coulomb branch, which in turn implies that the anti-D5 
branes' world- volume scalar multiplets, transforming un- 
der a certain subgroup of SU{N + M), will be responsible 
for the Higgsing mechanism. 

For the non-supersymmetric theory, this is rather easy 
to demonstrate. All we require is that the Higgs field 
should only transform under a certain subgroup of the 
first SU {N + M) group. The Lagrangian is: 



(2.5) 



where i = 1,2 refers to each SU{N + M) copy in the 
product gauge group. D^^cpk = d^4>k - igiA'l^XTi)ki4>i 
with (7i, Ai^, and Ti being the gauge coupling, gauge 
field and generators of the first SU{N -I- M) group re- 
spectively. We can choose the matrix representation of 
the Ti generator properly to demand what subgroup of 
SU{N + M) we want. 

Now we suppose the potential V{<j)) is minimized at 
= Vi. Then a generator Tf is broken if {Ti)^^ Vj ^ 0. 
To see this, let's write: 



(\)i{x) = Ui -I- i?i(x) 



(2.6) 



where Hi(x) is a real scalar field. The covariant deriva- 
tive of is: 

^^0^ = dM^) - mA1^{T^hk [vk + Hk{x)] (2.7) 
and the kinetic term for becomes: 



1 



1 

y^iJii K^lJij 



(2.8) 



where — igiT^-Vj. It is obvious now that A"' will get 
massive if T°-jVj ^ and thus the gauge group is broken. 
In our case, we only want to break M of the generators. 
How this is done depends on the details of the potentials 
and the specific values of N and M . 

From the dual gravity, we expect to see M anti- 
D5 branes at r — oo so that the gauge theory is al- 
most conformal. As the radial coordinate decreases, the 



number of anti-D5 branes become M^, as given earlier. 
For r << ro we expect the number of anti-D5 branes 
to completely vanish so that the gauge group becomes 
SU{N + M) X SU{N) henceforth the cascading behavior 
begins. In figure 3 we have plotted the behavior of the 
function /(r) = MJM. 



FIG. 3: A plot of the function /(r) 



for ro = 5 in 



appropriate units, and for various choices of a. For r << ro, 
the function vanishes whereas it approaches unity for r > ro. 



The supersymmetric case follows the same line of ar- 
gument as above. The general Lagrangian now is: 



C 



1 



327ri 



TtVfWl 



h.c. 



(2.9) 



where ($, Wa 



\DDDaV) are the appropriate A/" = 1 



chiral and the vector multiplets with ((/)fc, AJJ) being the 
complex scalar and the vector fields in their respective 
multiplets, is a gauge invariant Kahler potential, W is 
a gauge invariant superpotential, r = is the com- 

plexified gauge coupling with i9-angle and tr f is the trace 
in the fundamental representation. The FI terms don't 
appear because they are forbidden by the non-Abelian 
gauge invariance. 

The scalar potential obtained by expanding the above 
Lagrangian in components is a sum of and terms. 
The F and the D terms are: 



dW 

d(j)k 



(2.10) 



where denotes the generator in the R representation. 
To preserve supersymmetry we must have Ffe = and 
Z?" = 0. Actually in the absence of FI terms whenever 
Ffc — has a solution, = always has a solution. 
So we assume we already have a solution that satisfies 
Ffe = 0. This solution can break the gauge symmetry as 
in the non-supersymmetric case. This can be seen in the 
following way. 

Write down the relevant kinetic terms of the scalar 
component of the Higgs multiplet, as: 



implies that the anti five-brane degrees of freedom should only 
be seen at high energies, precisely in the way we predicted in Q! 



\(d^ + igAlT'^)<\>X + (2.11) 
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which is exactly the same as in the non-supersymmetric 
case. If the _F-term sokition leads to M generators T" 
such that T^-^Vj ^ 0, then the gauge group is broken from 
SU{N + M) to SU{N). How this happens again depends 
on the details of the iV, M values and the form of the 
superpotential. 



to the stress-energy tensor®. Thus, using the background 
above the deconfinement temperature given in [g, [gl, , 
the internal energy, in terms of the string coupling and 
the Newton's constant Gat, becomes: 

E,nt = / d'x^ Too - (3.13) 



3. PHASE TRANSITION AND OTHER 
APPLICATIONS 

Once we have the gauge theory description, it is time 
to extend our configuration to incorporate temperature. 
Two immediate scenarios present themselves: the confin- 
ing theory at low temperatures and the theory above the 
deconfinement temperature. The process of going from 
one to another in the gravity dual will appear as the con- 
finement to deconfinement phase transition in the large 
N thermal QCD. 

In [1, H, Q , the theory above the deconfinement tem- 
perature was studied in details. The high temperature 
phase was understood therein as the one coming from a 
black-hole with a horizon radius r/j where the tempera- 
ture was related to r^. The scenario at low temperatures 
were not discussed in details in ^, ^ Here, we will 
study these two phases and discuss their associate phase 
transition. More elaborations on this will be presented 
in m. 

Before actually computing the phase transition, let us 
discuss a couple of issues that may arise in equivalent 
scenarios dealing with large N thermal QCD. The first 
issue is the stability at high temperatures. Stability is 
guaranteed by a positive specific heat c„. A negative 
specific heat implies instability, which in fact turned out 
to be the case of many models that study large N thermal 
QCD without a UV completion [9|]. To assess the issue 
of stability, let us first define the specific heat in terms 
of the internal energy Eint, the BH factor g — 1 — rjjr'^ 
and the temperatute T in the following way: 



which gives the following value for the specific heat: 



fdE, 



int 



V dT 



T = 



(3.12) 



where we have introduced the AdS^ length scale L in 
anticipation of the AdS cap, and the internal energy is 
given by the integral of the zeroth component of the 
stress tensor. To calculate the heat capacity, we have 
to know how much energy is encoded in the geometry. 
As r — >■ oo, the space-time is approximately AdS^ x T^'^, 
where T^-^ is the internal space. The internal energy of 
asymptotically AdS ^-, sp ace-time can be easily calculated 
using results from [l3|. The total stress-energy tensor 
Tij is composed of stress-energy from the medium and 
the quarks. Quarks can be seen as excitations of the D7 
branes. At the boundary, only the medium contributes 



47r6L® , 

gs'^GN 



(3.14) 



This means that the heat capacity is positive for positive 
temperatures, showing that the model is stable at high 
temperatures. 

The second issue is slightly tangential to our interest 
but is nevertheless important enough that we clarify the 
scenario here. It was proposed recently in an interesting 
work [l6l | that the confinement to deconfinement phase 
transition in the type HA Sakai-Sugimoto model [13 does 
not proceed via the standard transition of a solitonic D4- 
brane to a black D4-brane, as proposed in J^, ilS], but 
via a Gregory-Laflamme transition (20l | from a solitonic 
D4-brane to a certain type IIB Euclideanized D3-brane 
configuration. In retrospect, this conclusion may not be 
too surprising because the black D4-brane elegantly de- 
picts the /lue-dimensional deconfined phase but fails to 
do so in the four-dimensional case once a certain energy 
scale is reached. Indeed in this model, there is no reason 
for integrating out the modes coming from the compact 
direction. Thus the Euclideanized D3-brane phase 
should be preferred at high temperatures. Unfortunately 
however, because the Mandal-Morita [16] picture above 
the deconfined phase is not a configuration of black D3- 
branes, the usual computations of transfer coefficients, 
that rely on the dynamics of black-holes in these spaces, 
cannot be performed so easily. 

This is exactly where our model may have some dis- 
tinct advantages. Since we are considering configurations 
of wrapped five-branes and anti five-branes on vanishing 
two-cycle, the subtlety of Kaluza-Klein (KK) reduction 
will not appear, and we should be able to go between 
solitonic D3 and black D3-branes. This would then be 
the confinement to deconfinement phase transition for 
our case, which is of course the Hawking-Page [22| tran- 
sition. In the following we will first take a brief detour to 
explain the Sakai-Sugimoto limit of our model, before go- 
ing into the discussion of phase transition in our set-up. 
More details will appear in [ll|, [2^ . 



* The contribution from Too a-iid warp factor /i as r — > 00 is of 
zeroth order. When the AdS geometry is deformed, the stress- 
energy receives higher order corrections. 
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3.1. The Sakai-Sugimoto limit 

From the above discussion, an interesting question at 
this stage would be to compare our type IIA dual pic- 
ture with the Sakai-Sugimoto model 17] ■ For simplicity, 
let us only consider the far IR picture where we have 
D5-branes wrapped on the vanishing two-cycle of the 
conifold. The vanishing cycle could be parametrised by 
(01, and the other two-cycle is along (02,</'2)- The 
[/(I) fibration of the conifold is and the radial direc- 
tion is r. The D5-branes have a spacetime stretch along 
the usual 2:°'^'^'^ directions. T-dualising along the ip di- 
rection gives us D4-branes stretched between two NS5- 
branes along the V' circle [21]. Thus the ij} coordinate 
is like the x'^ coordinate of the Sakai-Sugimoto model. 
The difference now is that the D4-branes are stretched 
only along a fraction of the if) circle and between two 
orthogonal NS5-branes. The D7 anti-D7-branes become 
D8 anti-D8-branes along r) and Pj^ x Pj^ 

just like the Sakai-Sugimoto case, but with an instanton 
configuration on the two spheres that breaks the super- 
symmetry. Of course, one might now worry that, since 
we made non-contractible, the usual issue raised in (T6| 
should appear for our T-dual model too. However, note 
that the distance between the two NS5-branes could be 
made arbitrarily small^ (without changing the size of the 
■0 circle), so the issue raised in 16^ may appear only at 
very high temperatures! Thus even at arbitrarily high 
temperature, if we tune the distance between the two 
NS5-branes appropriately so that the V' modes are of 
very high energies, we might still be able to study the 
deconfined limit using the black D4-branes. Further de- 
tails and explicit computations on this construction will 
be reported in 111). 



tifying it with the free energy of the gauge theory. 
We start with the on-shell type IIB supergravity action 
with appropriate Gibbons-Hawking boundary terms and 
counter terms; 



S — l3Efrcc — Sub + Sgh + Sc. 



(3.15) 



where Ehcc is the free energy. Sub is the ten dimen- 
sional type IIB Euclidean supergravity action including 
localized sources ^]^], Sgh is the Gibbons-Hawking 
surface term [2^ and ^counter is the counter term nec- 
essary to renormalize the action [13] [Bl] l23l] . Just like 
the case for AdS gravity discussed by Hawking and Page 
[2^ and subsequently by Witten Ts], the above action 
gives rise to both extremal and non-extremal metric and 
both geometries can incorporate non-zero temperature of 
the dual gauge theory in the following way: Wick rotate 
t — )■ ZT, T e (0, /3) and identify temperature T as T = 1//3. 
At a fixed temperature of the gauge theory, we have 
two geometries — extremal and non-extremal — and the 
geometry with smaller on-shell action will be preferred. 
The free energy of the gauge theory will then be given by 
the free energy of the geometry obtained through (|3.15p . 
Denoting the on-shell value of the action for the extremal 
geometry with Si and the non-extremal geometry with 
52, we compute the action difference in the absence of D7 
branes and localized sources, i.e. Nf = and the axio- 
dilaton r is a constant (i.e without fundamental matter), 
as ^23.]: 



AS 



S2 — Si 



lim 



00 



-log 
32 ^ 



128 



(3.16) 



3.2. Phase transition 

Phase transitions of SU{N) gauge theory can be re- 
alized by spontaneous breaking of the center symmetry 
Zjv- In the confined phase, Z^v symmetry is preserved 
and its associated order parameter, a temporal Wilson 
loop, is zero (i.e (W) =0). In the deconfined phase, 
Zjv symmetry is spontaneously broken with (W) 7^ 0. In 
@, we computed (W) using the gravity description and 
showed that OKS-BH geometry with large black holes 
give (W) 7^ while the OKS geometry without black 
holes give (W) = 0. This indicates that extremal geom- 
etry is dual to confined phase while non-extremal geom- 
etry corresponds to deconfined phase. 

Here we will obtain the critical temperature for con- 
finement / deconfinement transition by computing the free 
enegy of extremal and non-extremal geometries and iden- 



® This depends on the choice of the Bns field on the vanishing 
cycle [2l| . 



where Vg is the volume of R^xT ■ , T ■ being the base of 



the conifold with approximate radius L = (gsN) 



1/4, 



'a' 



N, M are number of D3 and D5 branes, 7^ is the bound- 
ary value of r, and rh is the black hole horizon radius. 
Here d > is a constant independent of A^, M, gs and de- 
pends on the boundary values of derivatives of the metric 
[23j . In obtaining (|3.16l) . we have only kept terms up to 
linear order in ggM^ /N which is valid for N ^ gsM"^ 
and the exact form of S'counter, -Sgh is presented in [23| . 
The critical temperature is obtained by evaluating the 
critical horizon rf^ for which AS'(r^) ~ and the result 
is [23.]: 



O 



V ^ TT exp (f ) (g,A^)i/4V^ 

(3.17) 

where we have used the scaling TZ = L = (ggNY^'^^/a' — 
00. For T > Tc, AS* < 0, i.e the black hole geome- 
try has lower free energy and thus preferred, while for 
T < Tc, AS* > 0, i.e the extremal geometry is pre- 
ferred. For extremal geometry, one readily gets an en- 
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tropy s — — ^ g^*^" = 0, while for the black hole geometry: complete geometry. 



(3.18) 



at lowest order in ggM'^/N and 6 > is a constant 
independent of N,M,gs. Observe that when M = 0, 
AS — 0,Vr/i — i.e extremal and non-extremal action is 
equivalent for all temperatures of the boundary gauge 
theory. This is consistent with the field theory picture 
because the M = limit gives an AdSs x T^'^ geom- 
etry which describes a conformal theory. A conformal 
theory on x with circumference /3 for has no 
phase transition since the value of /? can be scaled away 
by conformal invariance (isj — i.e the vacuum phase is 
equivalent to the thermal phase. 

Even with M ^ 0, when we do not have any D7- 
branes — i.e we do not have any matter in the fundamen- 
tal representation^" — the confinement to deconfinement 
phase transition for the gauge theory mimics the first or- 
der transition in pure glue theory and is described by a 
Hawking-Page transition in the dual geometry. 

Observe that in deriving p.l6p , we defined the bound- 
ary r = TZ ~> oo, but did not explicitly add a UV geome- 
try. By adding counter terms iScountcr to the on-shell ac- 
tion, we subtracted the terms in Sub + Sgh that diverge 
at the boundary r = TZ, which is effectively choosing a 
particular UV completion. Explicitly adding an AdS UV 
cap would require taking account of the localized sources 
in the bulk in addition to the fluxes, and the exact on- 
shell action for a UV complete geometry is not known. 
However, the UV completion resulting from our regu- 
larization already gives us a first order phase transition 
with an exact result for the critical temperature and thus 
is already insightful. Furthermore, since confinement is 
an IR phenomenon, the critical temperature may not be 
extremely sensitive to the details of the UV completion 
and thus the Tc in (|3.17l) can even be relevant for the UV 

The field theory has bi- fundamental fields A;, Bj and in the far 
IR can be equivalently described by pure glue SU(M) theory. If 
Tc is very small, the confined phase consists of glue balls and the 



4. CONCLUSION AND DISCUSSIONS 

In this paper we have managed to tie up some of the 
loose ends of our earlier works 0, H, Q related to the 
gauge theory description of the UV complete geometry 
predicted in the gravity side. The RG flow from UV to IR 
at zero temperature shows how the conformal behavior 
in the far UV ties up with the conflning dynamics in the 
far IR. The intermediate-energy physics is more involved 
and will be elucidated in our upcoming work flTl where 
we will also discuss how to evaluate the spectrum of the 
theory. As an interesting outcome of the UV completion, 
we could see how the stability of our background could 
be justified. Furthermore phase transition and related 
IR issues appear naturally in our set-up. If we ignore the 
fiavor branes, our gravity description gives us a first-order 
phase transition. Further details on this will appear in 
|23| . In the presence of the flavor branes, the physics 
is slightly more involved and will be discussed in [ll|. 
We have also managed to compare our model with some 
of the other models that study large N thermal QCD 
and showed how certain calculations may become more 
tractable in our set-up. Whether this is true for most of 
the other details of large N thermal QCD remains to be 
seen. 
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deconfined phase consists of free gluons of SU{M). If Tc is large, 
the deconfined phase is best described by Ai,Bj fields. 
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